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Chandra observation of the fast X-ray transient 
ICR J17544-2619: evidence for a neutron star? 
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Abstract. IGR J17544-2619 belongs to a distinct group of at least seven fast X-ray transients that cannot readily 
be associated with nearby flare stars or pre-main sequence stars and most probably are X-ray binaries with wind 
accretion. Sofar, the nature of the accretor has been determined in only one case (SAX J1819.3-2525/V4641 Sgr). 
We carried out a 20 ks Chandra ACIS-S observation of IGR J17544-2619 which shows the source in quiescence 
going into outburst. The Chandra position confirms the previous tentative identification of the optical counterpart, 
a blue 09Ib supergiant at 3 to 4 kpc (Pellizza, Chaty & Negueruela, in prep.). This is the first detection of a fast 
X-ray transient in quiescence. The quiescent spectrum is very soft. The photon index of 5.9±1.2 (90% confidence 
error margin) is much softer than 6 quiescent black hole candidates that were observed with Chandra ACIS-S 
(Kong et al. 2002; Tomsick et al. 2003). Assuming that a significant fraction of the quiescent photons comes from 
the accretor and not the donor star, we infer that the accretor probably is a neutron star. A fit to the quiescent 
spectrum of the neutron star atmosphere model developed by Pavlov et al. (1992) and Zavlin et al. (1996) implies 
an unabsorbed quiescent 0.5-10 keV luminosity of (5.2 ± 1.3) x 10'^'^ erg s~^. We speculate on the nature of the 
brief outbursts. 
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1. Introduction 

It has recently become clear that there is a distinct group 
of fast X-ray transients with outburst durations of order 
a few hours and peak fluxes of order 10~^ erg cm~^s~^ 
(2-10 keV) that are not readily identifiable with other 
similarly active X-ray sources: magnetically active nearby 
stars (i.e., DY Dra, RS CVn or pre-main sequence stars) 
or superbursters (e.g., Kuulkcrs 2004; time profiles, lu- 
minosities and spectra are inconsistent). The high fluxes 
and the lack of nearby counterparts suggest high lumi- 
nosities which would indicate an X-ray binary origin. The 
first case, XTE J1739-302, was discovered by Smith et al. 
(1998). Accurate Chandra localization yielded an identifi- 
cation with a highly reddened probably supergiant O star 
(Smith et al. 2003a). The list has been steadily growing 
and now consists of at least seven objects: XTE J1739- 
302, SAX J1818.6-1703 (In 't Zand et al. 1998), V4641 
Sgr (In 't Zand et al. 2000; Orosz et al. 2001), XTE 
J1901+014 (Remillard & Smith 2002), AX J1749. 1-2733 
(Grebenev 2004b), IGR J16465-4507 (Lutovinov et al. 
2005; Neguerela et al. 2005) and IGR J17544-2619 which 
is the subject of this paper. For the first two systems there 
is suggestive evidence that the strong variability is due to 
a high-mass companion star feeding the compact object 



through a wind instead of an accretion disk, like in many 
high-mass X-ray binaries (Smith 2004). 

IGR J17544-2619 was first reported from INTEGRAL 
detections on Sep. 17, 2003, when it exhibited two flares 
(Sunyaev et al. 2003; Grebenev et al. 2003) and a sec- 
ond time on Mar. 8, 2004 (Grebenev et al. 2004a). Three 
XMM-Newton observations (Gonzalez-Riestra et al. 2003, 
2004) show that even outside large flares, the flux ap- 
pears to vary violently, with fluxes between an upper 
limit of 5x10"^** erg cm~^s~^ and 4x10""'^^ erg cm~^s~"'^ 
(0.5-10 keV). The first detection of IGR J17544-2619 was 
with the BeppoSAX Wide Field Cameras that detected 
hours-long flares in 1996, 1999 and 2000 with 0.5-10 keV 
peak fluxes of a few times 10"^ erg cm^^s^^ (In 't Zand 
et al. 2004). Semi- weekly flux measurements since 1999 
with the RXTE Proportional Counter Array (Swank & 
Markwardt 2001) point to a duty cycle of about 5% above 
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erg cm-2s-i(2-10 keV; In 't Zand et al. 2004). 



Rodriguez (2003) tentatively identified a counterpart 
in the 2MASS catalog, 2MASS J17542527-2619526, which 
is also in the USNO Bl.O catalog. Its brightness is B2 — 
13.9, R2 = 12.0 and = 8.02. Pellizza, Chaty & 
Negueruela (in prep.) identify the optical counterpart as 
a blue supergiant of spectral type 09Ib at 3 to 4 kpc dis- 
tance. The nature of the accretor is thus far undetermined. 
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Fig. 1. Map showing source detections around IGR 
J17544-2619. All error circles are for a confidence level 
in excess of 90% (the radius of the 68%-confident ROSAT 
error circle was multiplied by two) . The Chandra position 
is indicated by an asterisk inside the XMM-Newton error 
circle. Technically, since the Chandra and XMM-Newton 
positions are so close to the ROSAT error circle an asso- 
ciation to the ROSAT cannot be completely ruled out. 

Rodriguez-Riestra et al. (2004) noted the similarity with 
SAX J1819.3-2525/V4641 Sgr which is also a fast tran- 
sient that at one time revealed a giant flare of approxi- 
mately IQ-^ erg cm-^s"^ (Smith et al. 1999; Hjellming 
et al. 2000). Orosz et al. (2001; 2003) and Orosz (2002) 
determined a mass between 6.8 and 7.4 Mq (Ict confi- 
dence) from optical spectroscopy, revealing that it is a 
black hole (BH) candidate. This similarity prompts the 
question whether the accretor in IGR J17544-2619 is also 
a BH. 

In this paper we present a 20 ks Chandra observation 
of IGR J17544-2619. Thanks to its higher angular resolu- 
tion and, therefore, smaller background level, Chandra is 
able to probe deeper than XMM-Newton which improves 
the signal-to-noise ratio of the quiescent emission. Such a 
study may yield constraints on the nature of the compact 
object. If the luminosity would be above 10^^ erg cm^^s^^ 
and the spectrum a black body with a temperature of 
order a few tens to hundreds of eV, this would repre- 
sent strong evidence for a neutron star (NS; Rutledge et 
al. 1999; Garcia et al. 2001). Furthermore, the Chandra 
observations provide an improvement in the positional 
accuracy by an order of magnitude. Thus, the optical 
counterpart can be identified unambiguously. This is not 
the first Chandra observation of a fast X-ray transient. 
XTE J1739-302 was observed for 5 ks (Smith et al. 2003a, 
2003b). However, inspection of those data shows that the 
source was not in quiescence. 
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Fig. 2. Top: Chandra hght curve of IGR J17544-2619 from 
all photons extracted from 4'.'5 from the centroid. The 
photon rate is not corrected for pile up. The background 
count rate is 2.5 x 10"'* phot . Bottom: Light curve from 
trailed image. No background subtraction was applied. 
The photon rate during the first 6300 s, 1.55±0.01 s^^, 
is due to the background. 

2. Chandra observation 

Chandra observed IGR J17544-2619 on July 3rd, 2004 
(obsid 4550), with the ACIS-S CCD array (Garmire et 
al. 2001) in the focal plane and no grating. The CCD 
frame time is 3.2 s, the exposure time 19.06 ks. We ana- 
lyzed these data with CIAO version 3.2.1. The source is 
clearly detected. The point source image shows a heav- 
ily piled-up source with a characteristic hole at the cen- 
ter and readout trails. The average photon position ac- 
cording to the Mexican Hat wavelet method employed 
in CIAO tool wavdetect is a20oo.o = 17''54'"25!284, 
^2000.0 = -26°19'52'.'62 {l^^ ^ 3?23, 6" +0?33). The 
nominal uncertainty is 0'.'6. This position is 1'.'2 from the 
XMM-Newton position (see Fig. and 0'.'2 from the pre- 
viously identified optical counterpart (Rodriguez 2003). 

Figure [3 (top) shows the time profile at 100 s reso- 
lution of the rate of all 4191 observed photons within a 
4'.'5 accumulation radius from the centroid. Interestingly, 
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Table 1. Power law spectral fits to the 4 different states indicated in Fig. |21 Errors are for 90% confidence. 
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State A^H (10^^ cm"'') F Obs. flux^ Unabs. flux^ xi 

quiescence* 06^^ 5.9±1.2 1.4^o ^ x 10""* g.S^fg" x lO"'"* 38%+ 

rise* i i 8.4t^;i x 10-^=^ ^-^tti x 10"" j 

peak'' 1.36 ±0.22 0.73 ± 0.13 2.3li[j{!^ x 10"^ 2.69t^ x 10"'' 0.839 (144 dof) 

postflare^ T T '^-^^to U x 10"^° l-STtg-'^^ x 10"^° T 



^Fixed; ''Average flux in 0.5-10 keV (erg cm~^s~^); * These spectra were extracted from the point spread function with 
extraction radii of I'.'S and 5'.'0 for the quiescence and rise data, respectively; ^ These spectra were extracted from the readout 
trails; "'"This is the goodness of fit as determined through 10^ Monte Carlo simulations. It is the percentage of simulations with 
a goodness of fit parameter better than measured, based on the best fit model. 



the source shows two markedly different states. During 
the first 6300 s the source is very faint with a mere 26 
detected photons. Subsequently, it rises in two steps over 
3700 s to a level that is approximately 250 times brighter. 
Thereafter it decays down to a level 4 times fainter at 
which it remains for the 5000 s remainder of the obser- 
vation. To assess the true dynamic range (corrected for 
pile up), we made a light curve of all 2934 photons in the 
readout trails that are at least 25" from the centroid, em- 
ploying CIAO tool acisreadcorr. Such counting rates are 
not subject to pile up for photon rates below a few thou- 
sand s"'^, because in the readout mode the CCD pixels 
have an effective readout time of only 40 /iS. The result is 
provided in the bottom panel of Fig. |21 The peak rate is 
214.5 ± 2.1 s-i (background subtracted). This is 6.5x10^^ 
times larger than during quiescence. 

We used the 40 fis resolved trail data to search for a 
ms pulsar signal by constructing Fourier spectra for every 
CCD frame with a net exposure time of 41 ms, and aver- 
aging the spectra. We followed this procedure for all data 
as well as for only those pertaining to the peak and post- 
flare intervals. No periodic signal was found. The 3a upper 
limit on the amplitude of a sinusoidal signal is ~10%. 

We extracted spectra for four characteristic flux levels 
which are indicated by dashed lines in Fig. [5] (top) : quies- 
cence, rise, peak and tail. The first two spectra were taken 
from within circles centered on the source position with 
different accumulation radii. The extraction radius for the 
quiescent spectrum was chosen to be 1'.'5, which is small in 
order to minimize the background (0.7 background pho- 
tons are expected in this region) . For a soft quiescent spec- 
trum (see below) the encircled energy is still 90 to 95% for 
such a radius. The spectra for the latter two time inter- 
vals, where pile up is heavy, were taken from the readout 
trails. Quick inspection of the spectra reveals that there 
is an obvious change in the hardness of the spectrum be- 
tween quiescence and thereafter. All photons of the former 
spectrum are detected between 1.0 and 2.2 kcV while for 
the rising part and the peak part the fraction of photons 
above 2.2 keV is 72% and 74%, respectively. 

We modeled all spectra with an absorbed power law, 
using XSPEC vs. 11.3.1 (Arnaud et al. 1996). For the qui- 
escent spectrum we applied the Cash statistic (Cash 1979) 
because of the low number of photons; for the remaining 
spectra the chi-squared statistic was used, after rebinning 



energy channels so that each channel contains at least 
15 photons. The three non-quiescence spectra were fit- 
ted simultaneously, employing a single free hydrogen col- 
umn density, which parametrizes the low-energy absorp- 
tion (Morrison & McCammon 1983), and a single photon 
index T. The results are provided in Tabled The fitted 
value for A^h is, within errors, identical to the interstel- 
lar value of (1.4 ± 0.2) x 10^^ H-atoms cm-^ (Dickey & 
Lockman 1990). When leaving free A'h or F, the fit did 
not improve (i.e., xt did not decrease). 

The quiescent spectrum could also be satisfactorily 
fitted with a NS hydrogen atmosphere model following 
Pavlov et al. (1992) and ZavHn (1996), with N^i fixed to 
1.36 X 10^^ H-atoms cm~^. The fit was similar in quality: 
37% of lO'' Monte Carlo simulations based on the best-fit 
model had a smaller C statistic than observed. We find for 
a NS at 3.5 kpc with M = 1.4 Mq, R =10 km and B=0 G 
a temperature of 90 ± 10 eV. The unabsorbed 0.5-10 keV 
flux is (1.9 ± 0.4) X 10-^^ erg cm^^s-i. 

3. Discussion 

The observation was fortunately timed catching the source 
both in quiescence and outburst. The outburst phase 
provides an unambiguous identification and the quies- 
cence phase a constraint for the nature of the accretor 
as follows. Kong et al. (2002) reports Chandra ACIS- 
S spectra of four quiescent BH systems and find that 
these can all be best modeled with an (absorbed) power 
law with F between l.lOto^l and 2.28tgj^ (90% confi- 
dence level error margins). Other simple models also fit 
the data, except the black body model. The power law 
model is consistent with an advection-dominated accre- 
tion flow according to Kong et al. Similar results were 
obtained on other BH candidates (Tomsick et al. 2003). 
The quiescent spectrum of ICR J17544-2619 does not ad- 
here to this rule. That suggests that ICR J17544-2619 is 
not a BH candidate, but a NS. The quiescent flux cor- 
rected for absorption and extrapolated to 0.01-10 keV is 
(3.5 ± 0.9) X 10~i^ erg cm^^s-i, based on the NS model 
spectrum. This translates to a luminosity at 3.5 kpc of 
(5.2 ±1.3) X lO'^^ erg s~^, a reasonable value for a NS. 

However, OB stars also emit X-rays with similarly soft 
spectra. In a list of optically bright OB-type stars ob- 
served with ROSAT (Berghhofer et al. 1997a and 1997b) 



4 



J.J.M. in 't Zand: Chandra observation of ICR J17544-2619 



there are 14 09I/II supergiants. Six were detected with 
luminosities between 1.7 and 7.1 x 10^^ erg s^^ and 
eight were not detected with upper Umits from 0.4 to 
1.7 X 10^^ erg s~^. Thus, a considerable portion of the 
observed luminosity may be explained by radiation from 
the donor star. Given the statistics, though, it is more 
likely that most comes from the accretor. 

We speculate on the origin of the outbursts, that are so 
brief compared to typical transient X-ray binary outbursts 
of at least a week. Our study shows that this is probably 
unrelated to the nature of the compact object, given that 
IGR J17544-2619 probably contains a NS while V4641 Sgr 
contains a BH (Orosz et al. 2001). Rather, it must be 
related to the donor star. As already noted by In 't Zand 
et al. (2004) and Smith (2004), optical counterparts for 
three fast X-ray transients have been identified and all 
are of early spectral type (09Ib for IGR J17544-2619, 08 
for XTE J1739-302 and B9 for V4641 Sgr). Many early- 
type stars are suspected to have massive winds driven by 
resonance line scattering that are highly structured and 
variable (e.g., review by Owocki 1997; see also Lamers & 
Cassinelli 1999). A star that is particularly interesting in 
the present context is the BO main sequence star t Sco 
whose otherwise homogeneous wind is thought to contain 
at any time 10'^ clumps with masses of order iQ'^^-'^o g 
(Howk et al. 2000). A sizeable fraction of these may fall 
back to the star, giving rise to shock X-ray emission and 
redshifted spectral features, while others may flow with 
the wind. Thus, it is fairly well established that clumpy 
winds can exist in early-type stars and clump capture by 
a nearby compact object may explain the time scale of 
the X-ray flares. Surprisingly, the masses of the clumps 
in r Sco are a good match to the fluence of the transient 
outburst within an order of magnitude, assuming that all 
gravitational energy provided by the fall of the clump on 
a NS is liberated in the form of radiation. Modeling of 
this effect for a supergiant rather than main sequence star 
and for an earlier spectral class than for t Sco should be 
carried out to verify whether the outburst duration and 
recurrence time make sense in this scenario. 

Our improved positional accuracy (by a factor of 
6 in one axis) confirms the identification with 2MASS 
J17542527-2619526. An initial spectroscopic study of that 
star is being carried out by Pellizza et al. (in prep.). As a 
follow-up it is important to carry out time-resolved spec- 
troscopy of this star. Measurements of Doppler-shifted 
spectral lines may provide constraints on key parameters 
such as the orbital period and the mass of the compact 
object. The shape of the lines is instrumental in deter- 
mining the outburst mechanism: it may show signs of the 
dumpiness of the wind (e.g., Howk et al. 2000) or of strong 
variability in wind speed and/or density (e.g., Lindstr0m 
et al. 2005). Preferably, there should be simultaneous op- 
tical and X-ray observations. 
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